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The Effect of Strain on Thermodynamics of the Weakly 
First-Order Phase Transition 



Elastic matrix distortion around a growing inclusion of a new phase is an- 
alyzed and the associated contribution to the Gibbs free energy is consid- 
ered. The constant-composition transformation from the parent to product 
phase is considered within the frame of Landau theory of phase transitions. 
The volume misfit between the inclusion and matrix is assumed to originate 
from the transformation volume change coupled with the phenomenological 
order parameter. The minimization of free energy with respect to the vol- 
ume change and order parameter gives the dependence of Gibbs energy on 
the volume fraction of the product phase. The transformation proceeds in a 
finite temperature region with the equilibrium volume fraction dependent on 
temperature rather than at a fixed temperature as it would be expected for 
the first-order transition. The activation processes are shown to be irrelevant 
and the transformation kinetics is found to be fluctuationless. 
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Spontaneous strain is known to play an important role in many cases of the phase trans- 
formations in solids [|I],|2| . Due to the difference in crystalline structures and elastic properties 
of the co-existing phases the volume and shape misfits appear between the inclusion of the 
product phase and surrounding matrix of the parent one. An associated elastic strain both 
in the inclusion and in the matrix tends to relax this misfit by the cost of an additional elas- 
tic energy ||. The minimization of this elastic energy defines the shape and crystallographic 
orientation of the new phase particles Q . The elastic effects change qualitatively the phase 
diagram of heterophase systems and the phase rule is no longer valid since it does not take 
into account the elastic interaction between the phases which co-exist in a state of coherent 
equilibrium ||. The 'chemical' free energy alone can not provide the equilibrium fractions 
of the phases with different composition through well-known double-tangent construction 

II- 

The elastic effects are considered in the present paper for the case of the first-order 
solid-state phase transition without change in a chemical composition. Such a transition 
may take place in pure elements, fixed-composition compounds and in alloys with very slow 
diffusion kinetics. The stresses are shown to result in a finite temperature interval of the 
phase co-existence with the equilibrium phase fractions depending on temperature whereas 
this co-existence appears at the single temperature point in the absence of elastic effects. 
The transition kinetics is shown to be changed as well and the activation processes does not 
control the transition development. 

A usual approach |3j to elastic effects is to let the transformation proceed in the inclusion 
disconnected from the matrix, to let the 'zero-stress' transformation strain tensor to appear 
and than to bring the inclusion into an elastic contact with the surrounding matrix. The 
strain energy is then calculated using disconnected inclusion as an elastic reference state for 
the product phase. Though the new phase is sometimes considered to have different elastic 
moduli, calculations are mostly being done in the linear elasticity theory. However, when 
the transformation is sensitive to an applied external stress, a 'zero-stress' assumption is no 
longer valid and a non-linear dependence of the free energy of the new phase on the stress 
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level should be considered. 

The simplest possible way to analyze the transition strain is to use Landau theory of 
phase transitions |||9| taking into account the coupling of the phenomenological order pa- 
rameter with the strain tensor components. Such a coupling has been recently considered 
for the case of a proper ferroelastic (martensitic) phase transition [ID and both the effect 



of hydrostatic and uni-axial pressure have been studied. In the present paper I consider the 
elastic interaction in a heterophase system for a model case of the spherical inclusions of the 
product phase embedded into an isotropic elastic matrix of the parent phase. The quadratic 
coupling of the order parameter with the volume change that corresponds to a variety of 
systems (e.g. see is analysed. 

The remaining part of this paper is organized as follows. We begin with a linear elasticity 
analysis of the displacement field and the elastic energy associated with the volume change 
in an inclusion embedded into the isotropic matrix. Then the free energy difference between 
the infinite crystals of parent and product phases is considered within Landau theory and the 
coupling of the phenomenological order parameter with elastic strain is taken into account. 
The energy cost of the formation of a finite fraction of the new phase as inclusions inside the 
bulk crystal of the parent phase is than considered. After minimization of this free energy 
with respect to a volumetric strain, the resulting expression that depends on the volume 
fraction of the new phase is analyzed. Then the transformation kinetics is briefly discussed 
and the role of fluctuations is analyzed. 

Let us start with the spherical inclusion in an isotropic matrix and consider the bulk 
crystal of the parent phase as a reference state for the elastic energy calculation. A spherical 
shape of the new phase inclusion presupposes either that transformation strain does not have 
shear components or that such components disappear upon averaging over the ensemble of 
twinned domains of the new phase composing the inclusion. The latter situation takes place 
in many cases of martensitic phase transformations in metals. 

The radial displacement field u(r) appears due to the transformation and it has the 
following form determined by the elastic equilibrium conditions \12\ 



u(r) 



a±r for r < R 
a 2 r + ^ for R<r< R 



Here i? is the inclusion radius and R Q = Rj ^/v is the radius of the spherical domain of the 
parent phase attributed to the inclusion for v being the volume fraction of new phase. The 
first boundary condition states that the displacement field vanishes at R Q and the continuity 
of the displacement field across the interface implies the second boundary condition. Hence, 
one can get 

R s 

di — — , CI2 = — - and 02 
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where e = Tr(e) is a volumetric strain inside the inclusion. 

The elastic energy associated with the matrix strain and expressed per unit volume of 
the new phase is given by an expression 

where 7 = (4« )/(3-K" ) is the renormed ratio between the shear and bulk moduli of the 
parent phase. This energy corresponds to the stress that plays a role of the external pressure 
applied to the inclusion of the product phase. 

The elastic energy (per unit volume of the new phase) associated with volumetric strain 
inside the inclusion is 

AC« - ¥° e 2 

Total elastic energy that has to be added to the Gibbs energy difference between the 
phases is equal to 

= ^I±^. (1) 
2 1 - v K ' 

This expression is valid for a sufficiently small v when the inclusions of the product phase 
are well separated and the overlapping of associated spherical domains of the parent phase 
can be neglected. 



The free energy difference between bulk crystals of the parent and product phases which 
are related by the symmetry breaking phase transition can be expanded in the mean-field 
Landau theory || in powers of the 'order parameter' rj. The theory initially was developed for 
the second-order transition where rj is continuous at the transition point, however, weakly- 
discontinuous first-order transitions can be considered as well |S| . If the symmetry groups of 
both parent and product phases are known a priory then the scalar order parameter can be 
used and the Ginzburg-Landau expansion of the free energy difference per unit of volume 
has the form |§J§ 

A£ GL (T, V ) = I (T - T c )rf + f + ^> . (2) 

where T c is a critical temperature. Only the second-degree coefficient is supposed to depend 
on temperature and the equilibrium value of rj is determined by the minimization of AQ. 
Stability requires the highest order coefficient C to be positive and the third degree term 
B ^ implies the first order of the transition. 

Gibbs free energy (0) has two possible minima. One with rj = corresponds to a high- 
temperature undistorted phase stable for T > T c . For 

1 B 2 

T<T = T C + - — 
4 aC 

a low-symmetry phase exists with 



i 



B L /T -T\5 



The phase energies become equal at the temperature of the first-order transition T* = 
T c + where the order parameter jumps from the i] = to 

2 B 

overcoming the activation energy barrier 



1 B 4 
AQb ~ 324 C^- 
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To analyze an associated strain effect the coupling of 77 with the strain tensor e has to 
be considered and the Ginzburg-Landau expansion should include the elastic terms ([I]). In 
the simplest case the symmetry allows the volumetric strain eo = Tr(e) to be coupled with 
rj in the lowest order by the following term: 

Ag int (e ,r]) = D e r] 2 . (4) 

Adding this term to the elastic energy (|l|) and finding the minimum of the resulting expres- 
sion with respect to eo we obtain the dependence of the volumetric strain inside the inclusion 
on the phenomenological order parameter in the form 

= (5) 
K 1 + 7 



This leads to a renormed Ginzburg-Landau expansion of AQ in powers of rj ||10||: 



Ag{T, v ) = Ag GL (T, v ) + Ag el (e (r))) + Ag int (e (v) , v) 

^4 ' V K C I + 7J 1 ; 

We have considered the unstrained parent phase as a reference state for an energy cal- 
culations, thus, the energy cost for the formation of the new phase (||) should be multiplied 
by its volume fraction. Choosing the case of B < which implies a positive rj in the prod- 
uct phase, we can write the expansion of total free energy per unit system volume in the 
following form 

r 3 

Ag = —Ag 



with 



n = ~C, r = ^(T-T c ) and * - ^ 



C * ' B 2 y c/ T K C 

The dependence of Ag on £ and v is shown in Fig.|I] for some r, ip and 7. 
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The state of equilibrium is determined by the minimum of this free energy with respect 
to both ( and v 

dAg , dAQ , s 

^T^° and "^7 = ° ^ 
d 2 Ag d 2 Ag 

^T + W >0 (8b) 

> . (8c) 



d 2 Ag d 2 Ag ( d 2 Ag \ 2 

d 2 ( d 2 v \ dvd( ) 
Solving coupled equations (pal) we obtain 



+ (9) 



3 V 1 + 7- V, 
1+7-V 9(l+7-^) 2 r 
27 ~~ 2^ 4 (1 + 7) ^ ' 

The stability conditions (|8b|) and (|8c|) are satisfied for these ( and z/ if 



(10) 



T<T*= 2 -(l+ t ). (11) 

9 \ i+7-vy 1 ; 

In the ^ — > limit r* corresponds to the temperature of the first-order phase transition T* 
in the absence of stresses. An important difference is, however, that now the product phase 
does not correspond to the free energy minimum above this temperature and, therefore, is 
mechanically unstable. The energy dependence on ( is shown in Fig||for r < r*. There is an 
activation barrier separating the initial £ = state and the product phase that corresponds 
to minimum of free energy for ( ^ 0. However, there is no barrier in the free energy 
dependence on v shown in Fig.[3[ It means that there exists a transformation path in the 
((, v) phase space that does not involve any activation process. 

The classical phase rule allows the equilibrium co-existence of the parent and product 
phases of the same composition only at the fixed temperature T*. However, the phase rule is 
based on the assumption that there is no elastic interaction between the phases. The finite 
equilibrium volume fraction of the second phase v that corresponds to the minimal Ag in 
the present model appears in some temperature interval below r* and depends linearly on 
the dimensionless temperature r according to Eq.(^) which is illustrated in Fig.^. The 



system has to be cooled further down for the increase in u, i.e. for the transformation to 
proceed. 

The present model corresponds to the 'athermal' kind of martensitic transformation in 
metals when the transformation begins at some start temperature M s , but the parent phase 
still exists until the temperature goes down to Mf, a martensite finish point. The finish 
point Mf corresponds to v — 1 and the present model can not be used in this region because 
the inclusions of the new phase are supposed to be well separated. 

In the absence of stresses the transformation rate is determined by Arrenius factor for the 
activation energy barrier as well as by the nucleation and growth kinetics. The system can be 
supercooled and the transformation proceeds with a finite rate at any constant temperature 
below T* and can be completed at the same temperature in a finite time. In the present 
model the transformation does not involve any (thermal) activation process as there is no 
energy barrier. It means that fluctuations do not play an important role, supercooling is 
impossible and there is no Arrenius factor in the transformation rate which therefore should 
be very high. Such a 'non-activated' kinetics is a characteristic feature of the martensitic 
transformations. If the temperature interval between M s and Mf is narrow enough than for 
any real cooling rate transformation would be completed almost instantly. This is known to 
be the case for so-called 'isothermal' kind of martensitic transformations in metallic alloys. 

To conclude we have analyzed an effect of elastic stresses in the matrix on the thermody- 
namics and kinetics of the first-order phase transformation in an embedded inclusion. The 
coupling of strain with the transformation degrees of freedom has been considered in the 
frame of Landau theory of phase transitions. The transition takes place in some temperature 
interval rather then at a certain temperature and the transformation kinetics appears to be 
activationless. This model may describe martensitic transformations in some metallic alloys. 
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FIGURES 
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Dependence of AQ on the order parameter £ and the volume fraction of the new 
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FIG. 2. Free energy vs the phenomenological order parameter £. 
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FIG. 3. Dependence of AQ on the volume fraction of the new phase v. Equilibrium v for these 
r, ip and 7 is 0.117. 
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FIG. 4. Temperature dependence of the equilibrium volume fraction of the new phase v for two 
different values of the ratio of elastic moduli 7. 
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